Vacuum polarization around stars: nonlocal approximation 
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We compute the vacuum polarization associated with quantum massless fields around stars with 
spherical symmetry. The nonlocal contribution to the vacuum polarization is dominant in the weak 
field limit, and induces quantum corrections to the exterior metric that depend on the inner structure 
of the star. It also violates the null energy conditions. We argue that similar results also hold in the 
low energy limit of quantum gravity. Previous calculations of the vacuum polarization in spherically 
symmetric spacetimes, based on local approximations, are not adequate for newtonian stars. 

^> \ PACS numbers: 

o : 
o ■ 

(N . I. INTRODUCTION 

> : 

In quantum field theory in curved spaces, the mean value of the energy momentum tensor (T^^(x)) plays a main 
role. When evaluated around black holes or collapsing stars, it contains the information about Hawking radiation. 
As a source in the semiclassical Einstein equations, it is responsible for the quantum corrections to the metric. 
It is well known that (T Mt ,(x)) is a nonlocal object, in the sense that its value at a given point x may depend of the 
I \ values of the metric in the past null cone of x (it may also depend on the topology of the manifold) . This is clear for 

■ example in cosmological situations, where the number of particles created depend on the whole history of the scale 
^0 ' factor. It is less obvious for collapsing stars and black holes (Hawking radiation only depends on the surface gravity 

. of the black hole), but in principle the exact (T^^x)} should depend nonlocally on the metric. 

Given a spacetime metric, it is in general difficult to compute {T^ v (x)) exactly. It is even more difficult to solve 
self-consistcntly the semiclassical Einstein equation with (T^(x)) as a source, since this would involve the knowledge 

■ of the vacuum polarization for at least a family of metrics. For these reasons, it is of interest to develop approximation 
\ and numerical schemes. If the quantum field is very massive, the DeWitt-Schwingcr approximation is adequate to 

describe the vacuum polarization Q . However, if the mass is small or vanishes, the situation is more complex. There 
q-i arc many papers dealing with this problem, starting from those concerned with conformal fields in Schwarzschild 
I i spacetime 0i 0, 0, @- Analytic approximations have been proposed for Schwarzschild, Reissner Nordstrom, and 
^Jy general static spherically symmetric spacetimes, for both conformal and non conformal fields |||. The approach of 
• • , contains as particular cases the previous analytic approximations 0, EL H to (T M „(x)) for black holes. It has been 
shown that the analytic approximation reproduces the exact result for Schwarzschild spacetime with high accuracy 
k>( | in the case of massless scalar fields (however this is not the case for spinor fields, see 0). 

5_j ■ The different approximations to (T M i,(x)) depend of course on the quantum state of the field. The states usually 
considered are the Hartle-Hawking, Unruh and Boulware vacua ; 8j . The Hartle-Hawking state describes a black hole 
in equilibrium with its thermal radiation. The Unruh state is the one that best mimics the gravitational collapse of a 
star. The Boulware state, being ill defined on the horizon, is unphysical for black holes. However, it is widely believed 
to describe the quantum vacuum around a non-collapsed object like a star or a planet. For these reasons, analytic 
approximations in the Boulware state have been used to compute vacuum polarization around stars [?| and to check 
the validity of the energy conditions [lj} ■ 

In the approximations mentioned above, the energy momentum tensor is determined by local functions of the 
metric and its derivatives. However, as we already pointed out, on general grounds one expects (T^(x)) to be a 
nonlocal function. In particular, in the weak field approximation, a covariant perturbative approximation to (T^^x)} 
is explicitly nonlocal [illll2j|. In this paper we will compute (T Ml/ (x)) around static spherical stars using this nonlocal 
approximation. We will consider massless fields with arbitrary coupling to the curvature in the Boulware quantum 
state. We will compare the result with the analytic approximations, and show that the nonlocal part dominates. 
Moreover, we will show that from the nonlocal result one can derive the well known 1/r 3 - quantum correction to the 
Newtonian potential 0, . Conversely, we will see that the nonlocal approximation can be easily understood from 
this quantum correction and the superposition principle, valid in the weak field approximation. 

The paper is organized as follows. In the next section, as a warm-up, we discuss the local and nonlocal approxima- 
tions to (4> 2 ). We expand the nonlocal approximation in a multipolar expansion, which shows explicitly its dependence 
with the internal structure of the star. We also show that there are surface divergences on the boundary of the star 
unless the Ricci tensor is sufficiently smooth. In Section III we compute (T^(x)) in the weak field approximation for 
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an arbitrary, static, spherically symmetric star. We point out again the dependence on the inner structure of the star 
and the existence of surface divergences. Section IV contains some applications: first we show that (T^ v (x)) violates 
both the null energy and the average null energy conditions. Then we compute quantum corrections to the exterior 
metric, and find that they depend on the internal structure of the star. Section V contains a short discussion and 
conclusions. 

Throughout this paper we use units in which h = c = 1, while retaining G = lz ^ 1. Our sign conventions are 
in the nomenclature of [Tal 



II. NONLOCAL APPROXIMATION FOR (<j} 2 ) 

In this section we compute the quantity ((j) 2 }, for a massless scalar field with arbitrary coupling (£) to the curvature, 
in a weak background gravitational field. We assume the classical source that generates the field to be static and 
non-relativistic, i.e., its stress energy tensor takes the form T^x) = p(x)5"(5°. 

An expression for ((f) 2 ) can be obtained by taking the coincidence limit of the Fcynman propagator, 

(<p 2 (x)) = -Im ( lim G F {x,x')^, (1) 
and then renormalizing it. With this purpose we first find the Feynman Green function by solving the equation 

[n x +ZR(x)]G F (x,x') = ±=5{x-x'). (2) 

Since we are working in a weak gravitational field, we assume that both the metric and the Green function differ 
only slightly from their Minkowski space counterparts: g^ u = r}^ + h^ v and Gf{x,x') = G^p\x : x') + G p (x, x'). 
Expanding Eq. @ in these small quantities and transforming Fourier in x, it is straightforward to solve it to first 
order as 



G F (x,x>) = GP(x,x') + -^ J d 4 fc J d 4 fc'y ^x^^\iR{x)~h^di y - 2 , (3) 

where = h^ — ^rj^Trft), and now R(x) is the Ricci scalar calculated to first order in h^. After changing variables 
to p = k — k' and q = k + k' , the q- integrals in Eq. (|3J can be performed by standard dimensional regularization 
techniques. Integrating in d-dimcnsional g-space and taking the coincidence limit leads to 



1 r r r AHx—x) r 



ik' (x — x') 



j j ln a R 

87r 2 4-d 256tt 6 p 2 . 



(4) 



where p is an introduced arbitrary mass scale. The first term in the square bracket in Eq. 0} is the divergent part 
of the Hadamard propagator which must be removed from the result to obtain a finite -and renormalized- value for 
G F {x, x). The second term, on the other hand, contains two pieces: one local and arbitrary proportional to R{x) In p? 
and with no physical significance, and other non-local and with physical meaning to (4> 2 ). 

Once dropped the divergent part in Eq. 0}, the integration in t and p° for a static situation yields 



f 1 / 1\ -V 2 
128tt 5 V s 6j J -r-vt" - p? ~ 128tt 5 y 6J n /i 2 

where R(p) is the Fourier transform of R(x). This result is finite everywhere for any smooth and asymptotically flat 
R(x). Notice, however, that for the external region of the source, i.e. R(x) = 0, Eq. JS} can be taken to a more 
friendly form if the d 3 p integral is performed first in Eq. 0}, namely 

(02(X)) = ~32^ ~ I) J i x ^|3 (external region). (6) 



As an application of the results in Eqs. i|5l6[) we now compute the vacuum polarization for a spherical star of radius 



i?o and constant density po (we assume pqGR 2 <C 1 to endorse the weak field approximation). In this case the scalar 
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curvature takes the form R(x) = 8ttGQ(Rq — r) and Eq. (JSJ) yields 



/^M\ G »°(t lA/M^-r 2 )] r<R 



2tt V" 6 J 1 -^ + ln(™) r > Ro, 

where C is an arbitrary constant proportional to p 2 . Notice that this expression has a logarithmic divergence at 
the surface of the star due to the discontinuity of the density; this is because the Fourier transform of the Heaviside 
function falls only as p~ l when p — > oo. In effect, any star model with a non-continuous p at the surface will have a 
logarithmic divergence in (cj) 2 ) at the point of discontinuity. On the other hand, for a continuous density the Fourier 
transforms falls faster to zero and the divergence disappears. The surface divergence is similar to the divergences that 
appear when computing the Casimir effect for perfect conductors of arbitrary shape. See for instance [16j . 

It should be noted that the expression for (</> 2 ) outside the star in Eq. (JJJ, depends not only on the mass of the 
star but also on its radius; in fact, expanding this expression for r > Rq yields 

MG / 3R 2 



^ 2 M> = -^^ — ( v 1 + 5^ + -"J (r>i?o) ' (8) 

where M is the total mass of the star. This kind of dependence is general for any star model, as can be seen by 
replacing Eq. © by its multipolar expansion, 



V)) = -— U--)f ^- 

1 47T 2 V 6 / ^ r 2n+l ' 

^ ' 71=1 



(9) 



which is exact in the external region for any spherically symmetric distribution of mass. The multipolar coefficients 
in Eq. are defined as 

i-Ro 

M n = 4tt / dfp(r)r 2n , (10) 
Jo 

so that Mi = M and in general M n ~ MRq . Therefore, the vacuum polarization in the external region depends 
at the sub-dominant order on the internal structure of the star. 

In contrast to this internal structure dependence, we point out that the local approximation for (0 2 ) in static 
spherically symmetrical spaces previously developed [T^ . for a massless field in the Schwarzschild metric leads to 

M 2 G 2 

W 1O0 = 4ta¥(2MGf _ r) . (ii) 

which, in the weak field approximation, goes as M 2 G 2 /r 4 since we are assuming r 3> 2MG in the external region. A 
comparison between Eq. and Eq. @ for large r gives 

{4> 2 (r -+ oo)) nonloc / 1\ r 

(0 2 (r)) Ioc e^MG' 1 > 

and hence, except in the conformal case ^ = 1/6, the nonlocal approximation is dominant. Moreover, since the nonlocal 
expression in Eq. Q increases when r — > R outside the star, the nonlocal expression is dominant everywhere in the 
external region. 

We have therefore shown through an easy object, as (</> 2 ), how the quantum effects on a weak gravitational field 
can induce dominant nonlocal behaviors, providing a conceptual difference with previous works in the subject. Our 
next step is to use the light thrown by this calculation to obtain and study analogous results for a more complicated 
and important object, the stress-energy tensor. 

III. NONLOCAL APPROXIMATION FOR (T„, .,) 

A nonlocal formal expression for (T, w ) in the weak field limit, similar in character to our expression for (0 2 ) in 
Eqs. Q-©, has been found in [Til. Il2|. In the massless case it reads 




l\ 2 _l_ 
90 



x\n^H$\x) + ±ln^H$>{x)\ , (13) 
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where the tensors H^l arc the two independent higher-order tensors that appear in the Einstein equations when 
derived from an action including R 2 and R^R^ . In the weak field limit they reduce to 

+ 0(R 2 ) . 



ff£> = 2V„V„ J R- 



tg^UR + O (R 
g^UR - 2UR 



(14) 



The action of the nonlocal operator F(D) = In in Eq. H13|) has been described in detail in several papers (see for 
example [TTI]'). For time dependent situations it involves an integral in the past null cone of x, and therefore (T,,, ,,) i s 
nonlocal and causal. On the other hand, for time independent situations, it can be shown that F(D) = F(— V 2 ) |18j . 

From Eq. i|13[l . assuming a static situation and performing the relevant integrations, it is possible to derive a 
nonlocal expression analogous to ©, valid in the external region where the classical source vanishes. It is given by 



(7^0r)> 



1287T 3 



1 

90 



d 3~ *w W 



- |3 
x — X 



1 

30 



d 3 i 



xl 3 



(15) 



Inserting the definitions Eq. (|14() into this equation, after integration by parts we obtain 
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(Too) = 
(Toi) = 



i 



32tt 3 




1 



32tt 3 



f-5 



1 
90 



1 

180 



V 



d 3 i 



R{5;) 



{didj - <\X 2 ,) I <f.r 



R(x) 



(16) 



therefore (T M1/ ) can be written in terms of derivatives of (<^ 2 )- For a spherically symmetric situation, using the series 
expansion Eq. (J5J) we find 



90 



^n(2n+ 1) 



+ I80 



M„ 
„2n+3 



M„ 



E( 2n+1 ) 

n=l 

o 



(17) 



Note that the nonlocal contribution to (T^) does not vanish for any value of £. 

The components of (T^) fall off as MG/r 5 for r ;§> R<j. On the other hand, in the local approximation for the 
Boulware vacuum (see for instance 0,^3), they fall off as M 2 G 2 /r 6 . So the nonlocal approximation dominates for 
(T^) as well as for (</> 2 ), by a factor of order r/AIG. 

Inside the star, (T M1/ ) can be found by taking derivatives of the finite expression for (0 2 ) (Eq.JSJl) in the same 
fashion. As an example, the complete (l 1 ^) for a star of constant density is, for r < Rq, 



<n> = 


Gp Q 

TT 


(T r r ) = 


2Gp 




TT 


<n> = 


2Gp 

TT 



1 

90 



+ T80 



3i? 2 - r 2 
(i? 2 - r 2 ) 2 

1 



7? 2 



1 

180 



^0 

(R 2 -r 2 ) 2 ' 



(18) 
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while for r > Rq 



7T 



^ 6 J 90 



r(r 2 



2Gp 



2Gp 



i 

180 



Rl) 2 
Rq 



i( r 2 



1 

180 



i?g(3r 2 



2r 3 (r 2 



(19) 



This result is quadratically divergent on the surface of the star. To remove the logarithmic divergence in (<fi 2 ) it was 
enough to require R(x) to be a continuous function; as (T^) is constructed with second derivatives of (</> 2 ), to ensure 
that no divergence arises on the surface it is necessary to require the continuity of the second derivatives of R(x). 

The local approximation for (T) t „) inside a star of constant density, presented in [jj, gives an expression with an 
additional factor of MG/Rq. Therefore the nonlocal approximation is, for nonrelativistic stars, dominant over the 
local one both in the internal as well as the external region. 



IV. APPLICATIONS 
A. Energy conditions 



One important question to ask is whether (T^) calculated within semiclassical gravity satisfies the energy conditions 
or not. It is well known that quantum fields often do not satisfy the energy conditions |19|, at least not in the same way 
as classical fields do. This fact has raised considerable discussion regarding its implications for singularity theorems, 
black hole dynamics, existence of macroscopic traversable wormholes, creation of closed time-like curves, etc |20j . In 
other words, violation of the energy conditions imply that the semiclassical Einstein equations, with {T^ v ) as a source, 
could in principle admit solutions qualitatively different from classical solutions |'2l(. 

In this section we will discuss the validity of two energy conditions, the Null Energy Condition (NEC) and the 
Averaged Null Energy Condition (ANEC) |22|. for the massless scalar field (T^) found in last section. The discussion 
will be restricted to the exterior region, since inside the star all energy conditions are satisfied by the classical source. 
These conditions have been studied in Ref. |Toj using the local approximation for the Boulware vacuum in the 
conformal case. The conclusion was that neither of them holds in the exterior of a non-collapsed star. However, since 
the nonlocal contribution to (T^) dominates in the weak field limit, it is necessary to revise the calculation. 

The NEC states that T ilv K ll K v > for every null vector K^. For a static spherically symmetric situation this 
condition is equivalent to the following inequalities 

(n>>(n> (T t t )>(T° e ) , (20) 

to which we will refer as radial and tangential conditions respectively. Using the series expansions (|17f) we find that 



(T\) ~ (T r r ) 



— Y 

2tt 2 ^ 



Mn 
„2n+3 



(2n + l) x 



1 

90 



n -~2 



e-g) (n+l) 



(21) 



(r* ) - <n> = - 



G ^ M n 

2tt-2 2-^1 r 2n+3 

71=1 



(2n+ 1) x 



360 



(6n+l) + 



1 

360 



(22) 



It is easily seen that the tangential condition is violated for any £, which is enough to ensure NEC violation. The radial 
condition is also violated in the conformal case, but it holds in the minimal coupling case. Therefore the possibility 
remains that in this case the average of the NEC condition over a null geodesic may be positive. We will show now 
that this is not the case: the ANEC is also violated for arbitrary £. 

The ANEC H3 states that dp T^K^K V > for any null geodesic of tangent vector and affine parameter 
p. To prove that the ANEC is also violated, we consider a null geodesic external to the star with impact parameter 
b > Rq, defined by the following functions x^(p) 



t{p) = P r {p) = Vp 2 + ° 2 
6{p) = n/2 tp(p) — arctan(p/£>) 



(23) 
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The geodesic is taken to be a straight line, to keep the analysis consistently at first order. The integration of 
(T^v) K^K 1 ' over the whole trajectory (with = ^-) can be performed using the series expansion for (T^) in 
Ea.lfmi. The result is 

G ^ M n 



/CO S~1 ~^ 

dp (T uv ) K^K V = — - V 
y K " w 2880 b 2 ^ 



J (26) 2 ™ 

1 , 2(4n(n+l)-3)r(2n + 2) + r(2n + 4) 
(2n+1) I> + 3/2)r(n + 5/2) • (24) 

This quantity is clearly negative and independent of £. We conclude that the ANEC is violated for a massless scalar 
field in the exterior of a nonrelativistic star independently both of the coupling of the field and of the internal structure 
of the star. 



B. Quantum corrections to the metric 



In this section we shall solve the linearized semiclassical Einstein equations, which take into account the backreaction 
of the quantum field over the spacetime metric. The results will be found to agree with those obtained by |l3j but 
disagree with those found relying on the local approximation for (7)^) 0- 

Disregarding both cosmological constant and higher-order terms, the semiclassical Einstein equations read 



(25) 



We take the source term to be separated in a classical part, the ordinary density of the star, and a quantum part, 
which is the (T^ v ) we have calculated in the previous section. The gravitational field is similarly separated in a 
classical term (the background spacetime) and a correction that is produced by the quantum (T^) 



We express the correction to the gravitational field in isotropic form. In Cartesian coordinates 



h iq) - 



h\f=g(r)S ij 



(26) 



(27) 



It is a simple matter to solve the equations ()25|l. linearized in h^ u , using the (T M!/ ) given by Eq. Ijl5|) as a source. The 
result is 



f(r) = 



G 
Air 2 



1 



1 

45 



d 3 x 



R(x) 



x — X 



(28) 



9{r) 





1 

90 







fir) 



(29) 



Note that the correction h^J is proportional to (</> 2 ) (except in the conformal case £ = 1/6, in which (</> 2 ) vanishes 
but h/w does not). Eq. l(2*5|) is correct only outside the star, where (</> 2 ) is given by Eq. ©. Inside the star, h$ has 
the same proportionality with the finite (<f> 2 ) g iven in Eq. JSJ). 

In the exterior region we can replace Ea. (|28J) by its multipolar expansion 



f(r) = 

7T 



E; 

n=l 



Mn_ 

„2n+l 



(30) 



Here we see an explicit dependence of the external metric with the internal structure of the star, which can be thought 
as a "quantum violation" of Birkoff's theorem. This is possible because the space surrounding the star, although empty 
at the classical level, contains a quantum vacuum energy given by (T^). 
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As an example, the quantum correction to the Schwarzschild metric external to a star of constant density is predicted 
by our nonlocal approximation to be 



fir) 



1 



1 



MG 



1 + -^ 

5 r 2 



(31) 



with g(r) given by Eq. (|29|l . 

The leading term in Eqs. 1(30(1 and 1)31(1 for r 3> i?o is the quantum correction to the Newtonian potential of a point 
mass, since it contains no information about the internal structure. Therefore the complete Newtonian potential of a 
point mass when the effect of a quantum massless scalar field is taken into account reads 



$(r) 



MG 



, 1 V 1 
^6 



MG 2 



(32) 



This result agrees with previous ones [l3j . It also has the same form that the long-distance quantum corrections to 
$ due to gravitons, calculated in the low energy limit of quantum gravity |l4j . By contrast, the correction to the 
Schwarszchild metric found in using the local approximation for (T),„) is 



goo(r ->■ oo) = 1 



IMG M 2 G 2 



60nr 4 



(33) 



which does not agree with previous results. 

It is worth noting that our result satisfies consistently the principle of superposition; if the gravitational field of 
the star is calculated as a sum of infinitesimal contributions of the form given in Eq. (|32f) . the result takes us back 
to Eq. In fact, Eq. (|28|l is correct even if there is no spherical symmetry, as can be seen by noting that the 00 

component of the linearized Einstein equations reads (in the Lorentz gauge) 



Dh {q) - 



-WttG 



(Too)-i(T\>) 



(34) 



Using the expressions Eq. (|16|) for (T M „), it follows immediately that h q '(x) is given by Eq. (|28|l (if the point x is 
located outside the sources, i.e., if R(x) = 0). 

As an application in a non-spherically symmetric situation, we outline here the calculation of the long-distance 
quantum correction to the Kerr metric by considering a rotating source. The classical T^ v has now off-diagonal terms 



related to the angular velocity of the star, T oi (x) = pur cos#(— shi tpx + cos tpy)i 
in turn causes (T^) to have, outside the star, the following components 



(2) 

This produces a nonzero my which 



(Tor 



rV 2 



240?H 



o , p(x')r' cos Q'(— sin ip'x + cos (p'y)i 
a x , 



(35) 



in addition to those included in Eq. l|16f) . which are not modified by the rotation. Inserting (Toi) in the semiclassical 
Einstein equations, we find that the long-distance correction to the off-diagonal terms of the Kerr metric are of the 
form 



hoi = 



G 2 



107rr 5 



(J x r) 



O 



G 2 JRq 



(36) 



where J is the total angular momentum of the source. Note that the information about internal structure is again 
contained in the sub-dominant term, which always falls off two powers of r faster than the dominant one. Once again, 
similar results hold for the graviton quantum corrections |23j| . 



V. CONCLUSIONS 



In this paper we computed (4> 2 ) and (T^{x)) for a massless quantum field with arbitrary coupling to the curvature 
in the geometry of spherically symmetric Newtonian stars. We obtained expressions for both quantities as multipolar 
expansions, valid in the weak field limit, that show explicitly the nonlocal dependence of the quantum effects. As 
a byproduct, we have demonstrated that previous local approximations to the vacuum polarization in spherically 
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symmetric spacetimes do not apply to Newtonian stars. Indeed, the nonlocal part is the leading contribution that 
overwhelms the local one both inside and outside the surface of the star. The results are divergent on the surface if 
the star model is not taken as sufficiently smooth. 

We have shown that outside the star (T^ v {x)) violates the NEC and ANEC energy conditions. A local approximation 
for (T^ u (x)) was used in previous proofs of these results, which stand now on firmer grounds. Our results hold for fields 
with arbitrary coupling, while previously only conformal fields had been considered. We also computed the quantum 
corrections to the metric, and found a "quantum violation" to Birkoff's theorem: the external metric depends not 
only on the mass but on all the multipolar moments of the distribution within the star. From the same principles it is 
also easily shown that, due to quantum corrections, the gravitational field inside a spherical shell is slightly different 
from zero. 

We have pointed out that the nonlocal quantum correction to the metric outside the star is a consequence of the 
1/r 3 modification to the Newtonian potential and the superposition principle. This fact can be used to argue that 
the nonlocal modifications to the metric will be present in any quantum theory of gravity, since gravitons also induce 
l/r 3 -corrections to the Newtonian potential. These corrections are, of course, extremely small. 
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